The critical behavior of frustrated spin models with noncollinear order 
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We study the critical behavior of frustrated spin models with noncollinear order, including stacked 
triangular antiferromagnets and helimagnets. For this purpose we compute the field-theoretic expan- 
sions at fixed dimension to six loops and determine their large-order behavior. For the physically 
relevant cases of two and three components, we show the existence of a new stable fixed point 
that corresponds to the conjectured chiral universality class. This contradicts previous three-loop 
field-theoretical results but is in agreement with experiments. 
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The critical behavior of frustrated spin systems with 
noncollinear or canted order has been the object of in- 
tensive theoretical and experimental studies (see, e.g., 
Rcf. In spite of these efforts, the critical behav- 

ior of these systems is still unclear, field-theoretic (FT) 
renormalization-group (RG) methods, Monte Carlo sim- 
ulations and experiments obtaining different results. 

In physical magnets noncollinear order is due to frus- 
tration that may arise either because of the special ge- 
ometry of the lattice, or from the competition of different 
kinds of interactions. Typical examples of systems of the 
first type are three-dimensional stacked triangular anti- 
ferromagnets (STA), where magnetic ions are located at 
each site of a three-dimensional stacked triangular lattice. 
Examples are ABX3-type compounds, where A denotes 
elements such as Cs and Rb, B stands for magnetic ions 
such as Mn, Cu, Ni, and Co, and X for halogens as CI, 
Br, and I. They may be modeled by using short-ranged 
Hamiltonians for TV-component spin variables defined on 
a stacked triangular lattice as 
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where J < 0, the first sum is over nearest-neighbor 
pairs within triangular layers (xy planes), and the sec- 
ond one is over orthogonal interlayer nearest neighbors. 
In these spin systems the Hamiltonian is minimized by 
noncollinear configurations, showing a 120° spin struc- 
ture. Frustration is partially released by mutual spin 
canting, and the degeneracy of the ground-state is limited 
to global O(A) spin rotations and reflections. As a conse- 
quence, at criticality there is a breakdown of the symme- 
try from 0(N) in the high-temperature phase to 0(A— 2) 
in the low-temperature phase, implying a matrix-like or- 
der parameter. Frustration due to the competition of 
interactions may be realized in helimagnets where a mag- 
netic spiral is formed along a certain direction of the lat- 
tice (see, e.g., Ref. Q). The rare-earth metals Ho, Dy 
and Tb provide examples of such systems. 

The critical behavior of two- and three-component 



frustrated spin models with noncollinear order is con- 
troversial. Many experiments (see, e.g., Ref. H]) are con- 
sistent with a second-order phase transition belonging 
to a new (chiral) universality class. This is partially sup- 
ported by Monte Carlo simulations (see, e.g., Rcf. [jl| and 
references therein). On the other hand, three-loop per- 
turbative calculations at fixed dimension d — 3 and 
within the framework of the e-expansion || indicate a 
first-order transition, since no stable chiral fixed points 
are found for N = 2 and N = 3. These three-loop anal- 
yses show the presence of a stable chiral fixed point only 
for N > N c with N c > 3: N c = 3.91 § and N c = 3.39 [|. 

To explain these contradictory results it has been sug- 
gested that these systems undergo weak first-order tran- 
sitions, that effectively appear as second-order ones in 
numerical and experimental works. This hypothesis has 
been supported by studies based on approximate solu- 
tions of the Wilson RG equations jl , and by Monte Carlo 
investigations || of modified lattice spin systems which, 
according to general universality ideas, should belong to 
the same universality class of the Hamiltonian ([!]), and 
which show a first-order transition. 

For larger values of N, all theoretical approaches 
predict a second-order phase transition, but there are 
still substantial discrepancies between Monte Carlo and 
three-loop FT calculations (see the discussion of Ref. || 
for N = 6). 

All these considerations show that a satisfactory theo- 
retical understanding has not yet been reached. It is not 
clear whether experiments are observing first-order tran- 
sitions in disguise or field theory is unable to describe 
these rather complex systems. Of course, one may think 
that the observed disagreement is due to the shortness 
of the available series, thereby calling for an extension of 
the perturbative expansions to clarify the issue. 

FT studies of systems with noncollinear order are 
based on the O(A) xO(M) symmetric Hamiltonian HQ 
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TABLE I. Coefficients b[f of the six-loop expansion of 



where <f> a (1 < a < M) are M sets of TV-component 
vectors. We will consider the case M — 2, that, for 
vq > 0, describes frustrated systems with noncollinear 
ordering such as STA's. Negative values of vq correspond 
to simple ferromagnetic or antiferromagnetic ordering, 
and to magnets with sinusoidal spin structures . 

For N = 2, which is the case relevant for frustrated 
two-component spin models, an e-expansion analysis in- 
dicates the presence of four fixed points: the Gaussian 
one, an XY fixed point, an 0(4)-symmetric and a mixed 
fixed point. Using nonperturbative arguments ||, one 
can show that the XY fixed point is the only stable 
one H among them. However, the region relevant for 
frustrated models, v > 0, is outside the domain of at- 
traction of the XY fixed point, which would imply a first- 
order transition. However, it is still possible that other 
fixed points are present in the region vq > 0, although 
they are not predicted by the e-expansion. For N = 3, 
one may easily show the existence of an 0(6) fixed point 
for v — 0, which is expected to be unstable Q. Accord- 
ing to the three-loop analyses of Refs. |||| no other fixed 
points arc found for tA = 3, which would imply that the 
transition is of first order as well. 

In order to investigate the existence of new fixed 
points, we have considered the fixed-dimension perturba- 
tive approach, extending the three-loop series of Ref. ||] 
to six loops. As we shall see, the results of our six-loop 
analysis are somehow surprising, contradicting most of 
the earlier FT works. Indeed, the analysis of the longer 
series provides a rather robust evidence for the existence 
of a new stable fixed point in the XY and Heisenberg 
cases, with critical exponents that are in agreement with 
the experimental results. 

In the fixed-dimension FT approach one expands in 
powers of the quartic couplings and renormalizes the the- 
ory by introducing a set of zero-momentum conditions for 
the two-point and four-point correlation functions. All 
perturbative series are finally expressed in terms of the 
zero-momentum four-point renormalized couplings u and 
v normalized so that, at tree level, u ~ uq and v ~ vq. 
The fixed points of the theory are given by the common 
zeros of the (3- functions /3 u (u, v) and f3 v (u, v). In the case 
of a continuous transition, when m — » 0, the couplings 
u, v are driven toward an infrared-stable zero u*,v* of 
the j3- functions. On the other hand, the absence of sta- 
ble fixed points is usually considered as an indication of 
a (weak) first-order transition. 

In Tables || and |n] we present the six-loop expansion 
of the /3-functions for M = 2, associated respectively 
with the rescaled couplings u = 3u/(16ttR 2 n) and v = 
3u/(167r), where R K = 9/(8 + AT). Since FT perturbative 
expansions are asymptotic, the resummation of the series 
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-8(95 + 41TV)/2187 
400(TV - 1)/2187 
-118(TV - l)/729 
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0.27385517 + 0.15072806 TV + 0.0074016064 TV 
(TV - 1) (-0.22580775 - 0.018235606 TV) 
(TV - 1) (0.26899935 + 0.020014607 TV) 
(TV - 1) (-0.11448007 - 0.0052835205 TV) 

(TV - 1) (0.0089444441 - 0.0016624747 TV) 

-0.2792572 - 0.1836675 TV - 0.021838259 N' z 
+0.00018978314 N 3 

(N - 1) (0.31461221 + 0.055501872 N - 0.00082480232 N 2 ) 
(N - 1) (-0.44030468 - 0.073595414 JV + 0.0014028798 N^) 
(N - 1) (0.2585304 + 0.037648759 N - 0.001217462 N 2 ) 
(N - 1) (-0.059220508 - 0.0045473385 N + 0.00073177586 N 2 ) 
(N - 1) (0.0068557417 - 0.00025180575 N - 0.00025443043 N 2 ) 
0.35174477 + 0.26485003 N + 0.045288106 N 2 
+0.00043866975 N 3 + 0.000013883029 N 4 

5.1 (N - 1) (-0.50696692 - 0.12967024 N - 0.0014771485 7V 2 
-0.000084746086 N 3 ) 

4.2 (TV - 1) (0.80562212 + 0.20236453 N + 0.0017578027 Af 2 
+0.00018710339 N 3 ) 

3.3 (TV - 1) (-0.59063997 - 0.14164236 N - 0.00017117513 JV 2 
-0.00019649398 TV 3 ) 

2.4 (TV - 1) (0.21149268 + 0.045057077 N - 0.00090030169 N 2 
+0.000094681347 TV 3 ) 

1.5 (TV - 1) (-0.040302604 - 0.0073126947 TV + 0.00040413932 TV 2 
-1.0675626 x 10~ 6 TV 3 ) 

0,6 (TV - 1) (0.0021121351 + 0.00061386793 TV + 0.000014729631 TV 2 

-0.000012344978 TV 3 ) 

7\0 -0.5104989 - 0.4297050 TV - 0.09535750 TV 2 

-0.004001735 TV 3 + 0.00003226842 TV 4 + 1.410456 x 10 -6 TV 5 

6.1 (TV - 1) (0.8994695 + 0.3004550 TV + 0.01442862 TV 2 
-0.0001545778 TV 3 - 0.00001094158 TV 4 ) 

5.2 (TV - 1) (-1.571308 - 0.5310390 TV - 0.02432978 TV 2 
+0.0002518739 TV 3 + 0.00002911732 TV 4 ) 

4.3 (TV - 1) (1.354181 + 0.4660641 TV + 0.01959649 TV 2 
-0.0001820084 TV 3 - 0.00003790796 TV 4 ) 

3.4 (TV - 1) (-0.6451180 - 0.2271141 TV - 0.007968847 TV 2 
+0.00009326566 TV 3 + 0.00002730379 TV 4 ) 

2.5 (TV - 1) (0.1858305 + 0.0666827 TV + 0.001642878 TV 2 
-0.00006134290 TV 3 - 0.00001228467 TV 4 ) 

1.6 (TV - 1) (-0.02753284 - 0.009832742 TV - 0.00004985576 TV 2 
+0.00001827182 TV 3 +4.881350 x 10" 6 TV 4 ) 

0,7 (TV - 1) (0.002029855 + 0.0003424570 TV - 0.00005976447 TV 2 
+3.644114 x 10~ e TV 3 - 1.527023 x 10" 6 TV 4 ) 



is essential to obtain accurate estimates of the physical 
quantities. For this purpose we studied the large-order 
behavior of the expansion in u and v at fixed z = v/u. 
For z = v/u fixed and M = 2, the singularity of the Borcl 
transform closest to the origin, Ub, is given by 



for AR 2N > z > 0, (3) 
for z < 0, z > 4i?2TV, 
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where a = 0.14777422... and R K = 9/(8 + AT). Moreover, 
we find that for z > 2R 2 n the Borel transform has a 
singularity on the positive real axis, which however is 
not the closest one for z < 4R 2 n- Thus, for z > 2R 2 n 
the series is not Borel summable. 

In order to determine the fixed points we use the same 
method applied in Ref. Il(]] to the analysis of the RG 
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TABLE II. Coefficients of the six-loop expansion of 



| i i i | i i i | i i i | 

31 



1 1 1 1 1 1 1 1 1 1 1 1 1 
41 




1 2 3 1 2 3 

FIG. 1. Zeroes of the /9-functions for N = 2 in the (u, v) 
plane. Pluses (+) and crosses (x) correspond to zeros of 
f3u{u, v) and j3v(u,v) respectively. 
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FIG. 2. Zeroes of the /3-functions for N = 3 in the (u, v) 
plane. Pluses (+) and crosses (x) correspond to zeros of 
Pu(u,v) and fa[u,v) respectively. 
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0.64380517 + 0.11482552 JV - 0.0068647863 JV 2 
-0.65001599 - 0.058188994 JV + 0.014412109 JV 2 
0.25727415 + 0.0050681184 N - 0.0053156417 JV 2 
-0.026488854 + 0.0013330934 JV - 0.00053805835 JV 2 
-0.76706177 - 0.17810933 JV + 0.00016284548 N' z 
-0.00070068894 JV 3 

0.95092507 + 0.14711983 JV - 0.013730732 N 2 
+0.001856844 AT 3 

-0.55504891 - 0.065767452 N + 0.0088371199 N 2 
-0.0016981042 JV 3 

0.14198965 + 0.01319996 N - 0.0006025682 N 2 
+0.00067387985 N 3 

-0.014727547 + 0.0017404281 N + 0.00015290902 JV 2 

-0.00011202681 N 3 

1.096535 + 0.31582586 N + 0.0094338525 N 2 
-0.00049177077 N 3 - 0.000086193996 N 4 
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+0.00032481932 N 3 + 0.00027008406 N 4 
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functions of the cubic model. We resum the perturbative 
series by means of an appropriate conformal mapping 
p"T[ that takes into account the large-order behavior of 
the perturbative series at fixed z and turns the origi- 
nal series into a convergent sequence of approximations. 
To understand the systematic errors we vary two differ- 
ent parameters, b and a, in the analysis. We apply this 
method also for those values of z for which the series is 
not Borel summable. Although in this case the sequence 
of approximations is only asymptotic, it should provide 
reasonable estimates as long as z < 4i?2jv, since we are 
taking into account the leading large-order behavior. 

In Figs. ^ and || we report our results for the ze- 
ros of the /3-functions, obtained from the analysis of 
the Z-loop series, I = 3,4,5,6. For each /3-function we 
consider 18 different approximants with b = 3, 6, . . . , 18 
and a = 0, 2, 4 and we determine the lines in the (u, v) 
plane on which they vanish. Then, we divide the domain 
< u < 4 and < v < 5 into 40 2 rectangles, marking 
those in which at least two approximants of each (3 func- 
tion vanish. No fixed point is observed at three loops, 
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TABLE III. Critical exponents for N = 2 and N = 3. Our 
results are labelled by FT. Experimental and Monte Carlo 
results are reviewed, e.g., in Ref. Jl). See also Refs. |l2| . 
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consistently with Ref. As the number of loops in- 
creases, a new fixed point — quite stable with respect to 
the number of loops — clearly appears for 
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1.9(1), 
1.8(1), 



4.10(15), 
3.00(15), 
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for N 



(4) 
(5) 



We have been conservative in setting the error bars: all 
zeros of the approximants with 3 < b < 18 and < a < 4 
lie within the reported confidence interval. Notice that 
the fixed points belong to the region in which the series 
are not Borel summable, but still satisfy v* /u* < 4i?2jv- 
Therefore, we expect our resummations to be reliable, 
and the stability of the results with respect to the order 
of the series confirms it. 

We then compute the eigenvalues of the stability ma- 
trix. They vary significantly with the two parameters a 
and b and turn out to be complex in most of the cases. 
Nonetheless, the sign of the real part of the eigenval- 
ues is always positive showing the stability of the new 
fixed points. A reasonable estimate of the exponent uj 
is however impossible. We also compute the critical ex- 
ponents, by estimating the corresponding six-loops series 
at the fixed point, following Ref. JlOfl . The results are in 
substantial agreement with the experimental and Monte 



Carlo (MC) estimates, see Table [II 



Finally, we compare the six-loop results with the criti- 
cal exponents computed to 0(1/ N 2 ) in the framework of 
the large-A^ expansion. For example, 



_ ! 16 1 / 56 640 



N 2 +0 \ N 3 



(6) 



We find v = 0.858(4) for N = 16 and v = 0.936(2) for 
N = 32, which compare reasonably with the estimates 
that one obtains from Eq. (|^), i.e. v — 0.885 for N = 16 
and v = 0.946 for N = 32. 

In conclusion, the extension to six loops of the FT ex- 
pansions solves the apparent contradictions between field 
theory and experiments. We find that new stable chiral 



fixed points exist for two- and three-component systems. 
The estimated exponents are in substantial agreement 
with experiments, whose conclusions on the nature of the 
phase transitions are thus confirmed. However, we note 
that first-order transitions are still possible for systems 
that are outside the attraction domain of the chiral fixed 
point. In this case, the RG flow runs away to a first-order 
transition. 
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